COMPUTATION OF NEARLY EXACT 3D ELECTROSTATIC FIELD IN 

GAS IONIZATION DETECTORS 
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Abstract 



The three-dimensional electrostatic field configuration in 
gas ionization detectors has been simulated using an effi- 
cient and precise nearly exact boundary element method 
(NEBEM) solver set up to solve an integral equation of 
the first kind. This recently proposed formulation of BEM 
and the resulting solver use exact analytic integral of Green 
function to compute the electrostatic potential for charge 
distribution on flat surfaces satisfying Poisson's equation. 
As a result, extremely precise results have been obtained 
despite the use of relatively coarse discretization leading to 
successful validation against analytical results available for 
two-dimensional MWPCs. Significant three dimensional 
effects have been observed in the electrostatic field con- 
figuration and also on the force experienced by the anode 
wires of MWPCs. In order to illustrate the applicability of 
the NEBEM solver for detector geometries having multiple 
dielectrics and degenerate surfaces, it has been validated 
against available FEM and BEM numerical solutions for 
similar geometries. 

INTRODUCTION 

Electrostatic forces play a major role in determining gas 
detector performance. Hence, a thorough understanding of 
electrostatic properties of gas detectors is of critical impor- 
tance in the design and operation phases. Computation of 
electrostatic field is conventionally carried out using one of 
the following options: analytical approach, finite-element 
method (FEM), finite-difference approach and boundary 
element method (BEM). While the first of these possibil- 
ities offers the most precise estimation, it is known to be 
severely restricted in handling of complicated and realis- 
tic detector configurations. The FEM is the most popular 
being capable of producing reasonable results for almost 
arbitrary geometries |2|. However, for the present task, 
the FEM is not a suitable candidate because through this 
method it is possible to calculate potentials only at spe- 
cific nodal points in the detector volume. For non-nodal 
points, interpolation becomes a necessity reducing the ac- 
curacy significantly. Moreover, numerical differentiation 
for obtaining field gradients leads to unacceptable electric 
field values in regions where the gradients change rapidly, 
for example, near the anode wires of an MWPC. The BEM, 
despite being capable of yielding nominally exact results 
while working on a reduced dimensional space(unknowns 
only on surfaces rather than volumes), suffers from inaccu- 
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racies near the boundaries. They also necessitate relatively 
complicated mathematics. 

In this work, we propose the application of a recently 
formulated Nearly Exact BEM solver (NEBEM) 1 1 1 for 
the computation of electrostatic field in gas detectors. The 
BEM is a numerical implementation of boundary integral 
equation (BIE) based on Green's function method by dis- 
cretization of the boundary only. For electrostatic prob- 
lems, the BIE can be expressed as 



m 



f G{ 

Js 



r,r')p{r')dS' 



(1) 



where (j){r) represents the known voltage at r, p(f^) rep- 
resents the charge density at f^, G{f^ f) = l/47re|r — f^| 
and e is the permittivity of the medium. In order to de- 
velop a solver based on this approach, the charge carry- 
ing boundaries are segmented on which unknown uniform 
charge densities (p) are assumed to be distributed. The un- 
known p and the known (j) are related through the influence 
matrix A 
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where Aij of A represents the potential at the mid point 
of segment i due to a unit charge density distribution at the 
segment j. By solving (O, the unknown p can be obtained. 
Major approximations are made while computing the in- 
fluences of the singularities which are modeled by a sum of 
known basis functions with constant unknown coefficients. 
These approximations ultimately lead to the infamous nu- 
merical boundary layer due to which the solution at close 
proximity of the boundary elements is severely inaccurate. 



PRESENT APPROACH 

In the present approach, we have proposed to use the 
analytic expressions of potential and force field due to a 
uniform distribution of singularity on a flat rectangular sur- 
face in order to compute highly accurate influence coef- 
ficients for the A and for calculations subsequent to the 
solution of p. By adopting such an approach, it is possi- 
ble to improve up on the above-mentioned assumption of 
singularities concentrated at nodal points and move to uni- 
form charge density distributed on elements. In general, 
the potential (/) at a point (X, F, Z) in free space due to uni- 
form source distributed on a rectangular flat surface on the 
XZ plane and having corners at (xi, 0, zi) and (^2, 0, Z2) 



is known to be a multiple of 
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Closed form expressions for the above integration and also 
for the force vectors have been obtained and used as the 
foundation expressions of the NEB EM solver described in 
in . It should be noted that these new foundation expres- 
sions are valid throughout the physical domain including 
the close proximity of the singular surfaces. In the present 
work, we have extended the NEB EM solver to solve the 
electrostatic problem of gas detectors as follows. 

For gas detectors, it is also be useful to model the elec- 
trostatics of a wire. For this purpose, we have considered 
the wires to be thin in comparison to the distance at which 
the electrostatic properties are being estimated. Under this 
assumption, the potential ^ at a point (X, F, Z) due to a 
wire element along the Z axis of radius r and length L is 
as follows: 



, , JX^ -^Y^^iZ- hf -(Z-h), 

where r is the radius of the wire and h is the half-length of 
the wire element. Similar expressions for the force field 
components have also been obtained which are not pro- 
vided here due to the lack of space. In addition to the ex- 
pressions given in Ol, (|4|l and the companion force field 
expressions have been incorporated in the NEB EM solver 
to compute the electric potential and field in ionization de- 
tectors. 



can be easily solved on modest desktop computer. Please 
note that the solver solves a complete three-dimensional 
problem. We have executed our codes on a Pentium IV 
machine with 2GB RAM running Fedora Core 3 and it took 
approximately half an hour of user time to solve the most 
complicated problem. 

RESULTS 

In this section, we will concentrate mostly on the electric 
field, rather than the potential, to save space. In Fig[l] we 
have compared the midplane electric fields computed by 
the present approach and analytic electric field computed 
by (2! for a 5mm x 5mm x 50mm larocci tube. It can 
be seen that the two results match perfectly well. However, 
when the length of the tube is reduced to 5mm, the three- 
dimensionality of the device becomes relevant even at the 
mid-plane. Thus, in Fig|2|the difference between the two 
values have become apparent. In order to illustrate the ef- 



^ 5X5X50mm^, 20|im 



Computed! 
Garfield ■ 



05 0.1 0.15 0.2 



Figure 1 : Comparison of electric field at the midplane 



NUMERICAL IMPLEMENTATION 

In this work, we are going to present results related to 
larocci tube, MWPC and RPC. These devices have flat 
surfaces as their boundaries. The larocci tube and the 
MWPC have anode wires, in addition to the flat conduct- 
ing surfaces. We have considered larocci tubes and MW- 
PCs of various cross-sections (5mm x 5mm to 16mm x 
16mm), lengths (5mm to 160mm) and anode wire diam- 
eters (20/im to 100/im). The anode wires have been as- 
sumed to be held at lOOOF whereas the other surfaces are 
assumed to be grounded. The flat surfaces in all these de- 
tectors have been segmented in to 21 x 21 elements. The an- 
ode wires have been modeled as cylindrical polygons of 12 
sides. Along the axial direction, these cylinders have been 
segmented in to 21 sub-divisions. These are once again 
flat surfaces and, thus, have been modeled by the same ap- 
proach as discussed above. The anode wires have also been 
represented as thin wires for which @ and related expres- 
sions have performed as foundation expressions. The max- 
imum number of elements with a wire representation for 
an larocci tube has been 1785 while that for a MWPC with 
five wires have been 1869. Under the assumption of cylin- 
drical polygons, these have increased to 2016 and 3024. 
These numbers are quite modest and the resulting system 
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Figure 2: Comparison of electric field at the midplane 

feet of three-dimensionality, we have presented the relative 
deviation defined as ^rror(%) = ^^^-^g^^^^^/^"^^ x 100 
for various larocci detectors in Figs|3land|4l In these fig- 
ures, results for wires represented as thin wires, as well 
as cylindrical polygons have also been presented. All the 
electric fields are at the mid-plane of the 3D detector, the 
length of which has been varied from 50mm to 5mm. It 
can be seen that the deviation becomes apparent (close to 
1%) as the detector length becomes twice the side of the 
square detector. When the length is the same as that of 



the side, the error is as large as 10% near the cathode sur- 
faces. In the latter case, it is more than 2% near the an- 
ode wire, which is huge considering the magnitude of the 
electric field near the anode wire. Fig. |4| illustrates the 
same point more clearly and, in addition, shows the dif- 
ference between the two representations of wire. It may 
be noted that electric field computed by the wire represen- 
tation seems to be consistently higher than that obtained 
using the polygon representation. 



larocci Tube : 5X5X50 mm , 20]im 
= 20M.m (1 0|j,m from anode wire) 



**' y = 2.490 mm (10|j,m from cathode surface) 
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Figure 3: Deviation of electric field at the midplane 
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Figure 4: Close up of deviation of electric field at midplane 

Next, we present the normal electric field computed 
along the axial direction of the larocci tube in Fig 13 The 
reference straight line has been provided by using 1 2 1 and 
is the analytic solution for a 2D tube. It can be seen that 
the 2D value remains valid for almost 85% of the detector 
along the axial direction. However, in the remaining 15%, 
3D effects are quite prominent and hence, non-negligible. 
There is one more very important point to be noted. Be- 
cause of the nature of FEM, it almost produces oscillatory 
results in potential and field (less in the former and more 
in the latter) near the cathode surface and anode wires. But 
these are the locations where the results need to be most 
accurate! By using NEB EM, it has been possible to pro- 
duce perfectly smooth results without any hint of oscilla- 
tion. This precision can only be attributed to the exact 
nature of the foundation expressions of this solver. This 
remarkable feature of the present solutions should allow 
more realistic estimation of detector behavior in any situa- 
tion. 

Similar validation and comparisons have been carried 



Figure 5: Comparison of normal field along the axial direc- 
tion 

out for an MWPC having five anode wires with a wire pitch 
of 2.5mm leading to similar conclusions. There is one fea- 
ture which is considered to be important mainly for MW- 
PCs, namely, the electrostatic force that acts on the anode 
wires in an MWPC, especially, its positional variation. In 
Fig|6| we present the variation of the horizontal force act- 
ing on the different anode wires in a five-wire MWPC as 
we move along the length of the wire. It should be men- 
tioned here that the edge wires in the presented case are of 
100/im diameter, while those inside are of 20/im. It is seen 
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Figure 6: Variation of horizontal force along the axial di- 
rection 

that the horizontal force becomes quite considerable as we 
move towards the edge of the detector despite the use of 
larger diameter wire as guard wires. 

In the following FigQ we have presented the normal 
electric field contour at the mid-plane of the detector. The 
sharp increase in the magnitude of this field is quite evident 
near the wire locations. 

Several detector configurations in current use such as 
MSGC, RPC, GEM have multiple dielectric configuration. 
They also have extremely thin charged surfaces such as the 
graphite coatings of RPCs. Such surfaces may necessitate 
treatment of degenerate surfaces. As a result, in the follow- 
ing, we have tried to validate the NEB EM solver by com- 
puting the electrostatic properties of such configurations 
and comparing them with FEM and Dual BEM (DBEM) 
results given in |3| where the problem geometry is dis- 
cussed in detail. In the following Tables [T] |2| and |3j we 
present the comparison among the potentials as computed 
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Table 3: Comparison for R = 10 



Figure 7: Surface and contour plots of the normal electric 
field 



by the above three approaches. It is obvious from these re- 
sults that the agreement is excellent despite the use of only 
11x11 discretization of the surfaces. Please note that in 
the tables R denotes the ratio between the permittivity of 
two slabs placed on top of another (upper/lower) and the 
locations are expressed in microns. 



Location 


FEM 


DBEM 


NEBEM 


24.0,16.5 


0.514489 


0.52181 


0.5247903 


6.5,12.0 


0.230147 


0.23801 


0.2398346 


22.5,6.0 


0.3638855 


0.34638 


0.3451232 


4.0,3.5 


0.1108643 


0.10623 


0.1058357 
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Figure 8: Electric field in an RPC 



Table 1 : Comparison for R 


= 10 


Location 


FEM 


DBEM 


NEBEM 


18.0,3.0 


0.1723103 


0.17302 


0.1740844 


4.0,9.0 


0.2809692 


0.27448 


0.2807477 


25.0,16.0 


0.6000305 


0.59607 


0.5991884 


5.0,17.0 


0.679071 


0.67492 


0.6785017 



Table 2: Comparison for R = 


0.1 


Location 


FEM 


DBEM 


NEBEM 


18.0,3.0 
4.0,9.0 
25.0,16.0 
5.0,17.0 


0.01741943 
0.0281006 
0.4883313 
0.5929200 


0.017302 
0.027448 
0.480640 
0.589690 


0.0171752 
0.0286358 
0.4828946 
0.5926387 



Finally, in FigO we have presented the electric field con- 
figuration of an RPC computed using NEBEM. The RPC is 
assumed to have 2mm gas gap enclosed by two 2mm glass. 
On the two outer surfaces of the glass, graphite coating has 
been applied. The upper graphite coating has been raised to 
8kV, while the lower one is grounded. The RPC is 10mm 
wide, enclosed on the left and right sides by glass and has 
a similar spacer in the middle. The variation of the elec- 
tric field along the vertical direction (Y) on the midplane 
at various X locations have been presented. It can be seen 
that the normal electric field rises up to around 40kV/cm in 
the gas gap of the RPC. 

CONCLUSION 

Thus, it may be concluded that 

• Precise computation of three-dimensional surface 
charge density, potential and electrostatic field has 



been carried out for gas ionization detectors using the 
Nearly Exact BEM (NEBEM) solver. 

• The NEBEM yields precise results for a very wide 
range of realistic electrostatic configurations includ- 
ing multiple dielectric systems. 

As future plans we have the following aspects 

• Optimization of the NEBEM solver. 

• Multiphysics nature of NEBEM needs to be explored. 
After all, the foundation expressions represent new so- 
lutions to the Poisson equation which is one of the 
most important equations in physics ! 

• Applications to other areas will be explored 
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